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ABSTRACT 

We have measured dramatic flux and spectral variations through the 0.193 d orbit of the optical counterpart 
of the unidentified 7-ray source OFGL J2339.8 — 0530. This compact object companion is strongly heated, with 
T e f f varying from ~ 6900 K (superior conjunction) to < 3000 K at minimum. A combined fit to the light 
curve and radial velocity amplitudes implies Mi « O.O75M0, M2 « 1.4M© and inclination i m 57°. Thus 
this is a likely 'black widow' system with a E » 10 34 ~ 345 erg s pulsar driving companion mass loss. This 
wind, also suggested by the X-ray light curve, may prevent radio pulse detection. Our measurements constrain 
the pulsar's reflex motion, increasing the possibility of a pulse detection in the 7-ray signal. 
Subject headings: gamma rays: stars - pulsars: general 



1. INTRODUCTION 

The Large Area Telescope (LAT) on the Fermi satellite has 
revolutionized our understanding of the 7-ray sky, with par- 
ticular success in detecting spin-powered pulsars and blazars. 
After only three months of observation, a preliminary OFGL 
'bright source list' of 205 objects detect ed at > 10er signifi- 
cance was presented in lAbdo et al.l d2009al) . Remarkably, over 
95% of these source s now have been associated with lower 
energy counterparts (lAckerman et alJl201 ll) . Moreover, com- 
parison with the 7-ray properties of the identifications make 
it possible to classify many of the remaining sources; likely 
blazars show significant variability, while likely pulsars are 
steady, but show strongly curved spectra, with cut-offs above 
a few GeV. 

At present, there are two OFGL sources at |6| > 1° 
with pulsar-like properties and no identification. One, OFGL 
J2339.8— 0530 was a ~ 12cr detection with steady, hard spec- 
trum emission. It has been the subject of both radio pulsar 
cou nterpart searches an d 'blind' searches for 7-ray pulsations 
(eg. lAbdo et alj|2009bl) . To date no pulsed emission has been 
seen and the sourc e remains unidentified in the latest catalog 
(lAbdo et alJl20TTh . There it is listed as 2FGL J2339.6-0532, 
localized to I = 81.357 ± 0.03, b = -62.467 ± 0.02, with 
a variability index of 15.7 (indicating steady emission) and a 
flux of 3.0 ± 0.2 x 10" n erg/cm 2 /s, with a hard T = 1.96 
spectrum showing 6a evidence of a spectral cutoff. It thus 
remains one of the best unidentified pulsar candidates. In- 
deed, Kong et al. (201 lab), have drawn attention to this 
source, finding that there was an optically variable star co- 
incident with a CXO source in the LAT uncertainty region, 
arguing that this was likely a 'radio-quiet' millisecond pul- 
sar (MSP) and suggesting that it may be an LMXB in qui- 
escence, similar to the LAT-as sociated L MXB/MS P transient 
FIRS T J102347.6+003841 ^Archibald etaDl20Tol: IrarrTeTaT] 
2010). We report here on optical imaging and spectroscopy 
of this source which support the millisecond pulsar hypoth- 
esis, but instead show that this is a 'black-widow' type pul- 
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sar, evaporating a low mass companion, similar to the orig- 
inal of such systems, PSR B1957+20, but 2-3 x closer at 
d ~ 1 kpc. We suggest that strong outflow in this evaporat- 
ing system inhibits detection of radio pulsations. If the pulsed 
emission can be detected in the LAT photons, the resulting 
orbital information will make this a double-line spectroscopic 
binary and should allow an unusually precise determination 
of the neutron star mass, comparable to the important high 
Mpsr = 2.40 ± 0.12M W determination recently made for 
PSR B 1957+20 ( Ivan Kerkwiik. Breton & Kulkarnill20Til) . 

2. OBSERVATIONS 

Examination of plots of white light CCD photometry made 
on Oct. 31 and Nov. 1, 2010 with the lm Lulin Telescope 
(Taiwan) and on November 11, 2010 with the 0.81m Tenagra 
Telescope (Arizona) allowed us to infer an orbital period of 
^0.19d. The source appears in the SDSS DR8 data release 
with colors u=20.85, g=19.00, r=18.61, i=18.25, z=18.23, 
suggesting a G-type spectral class. Optical extinction at this 
high latitude is very small (Ay ~ 0.1 from the Schlegel et al. 
1998 dust maps). Based on the large (> 2.5mag) variation it 
seemed likely that this is a nearby, short-period black-widow 
type pulsar, with a strongly heated companion and thus a suit- 
able target for detailed optical spectroscopy. 

2.1. HET, SSO and WIYN Photometry 

The 9.2 m Hobby*Eberly Telescope (HET) has a spherical 
primary with a tracking corrector and can follow a source 
for ~ lhr/night during a transit. This is a small fraction of 
the estimated 4.6 h orbital p eriod. However the H ET is dy- 
namically queue scheduled dShetrone et al.l 120071) . and with 
an ephemeris one can obtain full orbital coverage. During 
HET Low Resolution Spectrograph (LRS) observations one 
obtains short direct images for target acquisition. We used 
these 3-15 s 'pre' images, augmented by 'post' images in sev- 
eral cases, to monitor the flux of J2339— 0533. These im- 
ages, through a Schott GG385 long-pass filter, approximate 
the 'white light' of the 2010 photometry. We measured sim- 
ple differential aperture magnitudes, calibrated to the SDSS 
r' magnitudes of nearby stars. These photometric data cover- 
ing Aug. 5 through Sept. 6 2011 confirmed the dramatic op- 
tical modulation of J2339— 0533. Conditions were variable, 
but detections were always of high S/N, even at the r' ~ 21 
minimum. 
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FIG. 1.— HET/SSO/WIYN measurements of J2339-0533. Upper panel: 
white light photometry (normalized to SDSS r), The curve gives the best-fit 
ELC model fit. Middle panel: WIYN SDSS band color photometry and ELC 
models. Lower panel: radial velocity measurements. Symbols indicates the 
MK class of the best-fit cross-correlation template. Curves show the model 
radial velocity for the Center-of-Light (solid) and Center-of-Mass (dashed) 
from the ELC fit. The dotted line shows 10 X the non-linear RVi variation 
10(Wi - 0.9/Cisin(i)). 

Since the target is very bright, we obtained additional pho- 
tometry with the 0.61m Cassegrain telescope at the Stan- 
ford Student Observatory (SSO) on Aug. 28, 2011 (MJD 
55801.287 - .399) and Oct. 20, 2011 (MJD 55855.140 - 
.381). Using 300 s unfiltered CCD exposures with an Apogee 
AP8 camera, we made differential flux measurements, again 
normalized to SDSS r' magnitudes. These data helped in ab- 
solute phasing and improving the binary period estimate. 

Finally, we obtained SDSS g'r'i' frames at the 3.6 m WIYN 
telescope using the MiniMo camera on Sept. 27, 2011 
(55832.290 - .424). These data, covering inferior conjunc- 
tion, were used to constrain the color variation through mini- 
mum (inferior conjunction IC). 

After converting exposure midpoint times to the solar sys- 
tem barycenter, we minimized residuals to obtain an accurate 
orbital period (Table 1, errors from bootstrap). The epoch 
(jointly fit with the radial velocity data, see below) defines su- 
perior conjunction (approximately maximum light). The pho- 
tometric data folded at this period are shown in figure 1 . 

2.2. Archival X-ray Light Curve 

We used this ephemeris to extract an X-ray light curve 
from the archival 2 Iks CXO ACIS-I exposure of the 0FGL 
J2339.8-0530 field taken on Dec 13, 2009 (MJD 55117.54, 
obs 11791, T. Cheung PI), covering 1.25 orbits. Figure 2 
shows the exposure-corrected light curve. The source is 
clearly modulated, with a relatively constant emission from 
phase -0.25 to 0.25, and a deep minimum at inferior conjunc- 
tion (phase 0.5). The source is, however, detected at mini- 
mum. 
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FlG. 2. — Two periods of the X-ray orbital light curve of J2339— 0533. 
Statistical errors on the bin exposure-corrected count rate are shown during 
the first period. A simple 1.5-8 keV/0. 5-1 keV hardness ratio is plotted (right 
scale) during the second period. 

A simple absorbed power-law fit to the phase-averaged 
data gives a hard index V = 1.09±of 3 and no significant 
absorption N H < 1.6 x 10 21 cm 2 with a (0.5-8keV) flux 
2 - 3 -o!3 x lCT^ergcm^s- 1 (all 90% errors). Lacking 
counts for a phase-resolved spectral analysis, we can extract 
a simple hardness ratio. Intriguingly, this is lowest at cj) = 0.5 
(IC), suggesting that the hard spectrum source dominates near 
SC, but is deeply absorbed or eclipsed revealing a softer com- 
ponent at IC. Deeper X-ray observations are clearly needed to 
probe the orbital spectral variations. 

2.3. Spectroscopy 

For the spectroscopic exposures we used the HET LRS, em- 
ploying grism 2, with a GG385 filter and a 1 .5 arcsec slit. This 
covered AA4283-7265A, at a resolution of 1.99 A/pixel, for an 
effective R ~ 1000 — 1200. All observations were taken with 
the slit at the parallactic angle and we limited exposures to 
600 s to minimize velocity smearing to < 35km/s. The HET 
pupil varies as the corrector tracks across the primary, so sig- 
nificant changes in the PSF and hence in the line profile of 
stellar sources occurs at the extrema of the tracks. Since we 
are seeking maximal radial velocity accuracy, velocities from 
observations very early or late in a transit may be compro- 
mised. For most HET visits, we were able to obtain 3-4- 600 s 
exposures; however a few exposures were taken far from the 
track center and were omitted from the orbital velocity fits. 
J2339-0533 lies relatively near the HET's DEC=-1 1° south- 
ern limit. It is thus possible to extend visit times by reset- 
ting the telescope azimuth during the exposure sequence. On 
Sept. 1, 2011, we obtained such a 'long track' re-setting the 
azimuth twice and allowing 9 x 600s exposure to be obtained 
(spanning ~ 2 hours or 0.42 of orbital phase) with only mod- 
est tracker excursion. Thus, under favorable circumstances, 
the HET can obtain extended integrations during a single visit 
at the cost of queue efficiency. In total 48 x 600 s exposure 
were obtained during one month of queue observations un- 
der variable conditions, equivalent to a full night of classical 
observation. 

Data reduction was p erformed with the IRAF package 
(Todv l[T986H \^ldes 1986) using standard techniques. Biases, 
dome flats and arc lamp exposures were obtained nightly. All 
data were bias subtracted, and flat fields were applied after 
removal of the lamp response. Arc exposures were used for 
wavelength calibration. Checks of the sky lines showed that 
these solutions were stable to ~ 0.3 A during a single visit 
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FIG. 3. — Averaged spectra for several phase bins during the orbit of J2339- 
0533. Green traces show three comparison spectra [F0,K2,M5]. The average 
of 0.8 < <f> < 0.1 at superior conjunction (maximum) shows an ~F-type 
continuum, with weaker Balmer emission and stronger metal lines. The av- 
erage of 0.4 < (f> < 0.6 at inferior conjunction (minimum) shows a M-type 
spectrum, dominated by the molecular complexes, but with some blue excess. 
Intermediate spectra represent phase ranges of 0.1 and are adequately traced 
by K-type spectra, trending toward M. 

and ~ 0.7A night-to-night; sky lines were monitored to sub- 
tract out these residual shifts, leading to an estimated stabil- 
ity of 0.1 A for the wav elength scale. Spectra were optimally 
extracted ( Valdesl [l992l) and visually cleaned of residual cos- 
mic rays. Spectrophotometri c calib ratio n was p e rform ed us- 
ing standard stars from lOkd {1990) and lBohlinl d2007l) . Oc- 
casionally appropriate standards were not measured during a 
given night; data from a neighboring queue night were then 
employed. The fluxing was not reliable for ~ 20A at the 
ends of the spectra and these were generally excluded from 
the final analysis. Spectra were corrected for atmospheric ex- 
tinction using the KPNO extinction tables. Telluric templates 
were generated from the standard star observations in each 
night, with separate templates for the oxygen and water line 
complexes. We corrected separately for the telluric absorp- 
tions of these two species. We found that most telluric fea- 
tures divided out well, with significant residuals only apparent 
in spectra with the highest S/N. 

The final cleaned, calibrated spectra were averaged to study 
the companion spectrum at various orbital phase bins (Figure 
3). Although the LRS does not cover much of the far blue 
used for traditional MK classification, there are clear con- 
tinuum and line dominance trends through the orbit. Com- 
parison with the Gray on-line atla£] and with sample spec- 
tra from the SPS S DR8 SEGUE Spectroscopic Parameter 
Pipeline (SSPPfl dLee. et al.ll2008h allows reasonably accu- 
rate temperature classification. The spectrum evolves from a 
~ 6900 K F3 type continuum at flux maximum to a ~ 2900 K 
M5 class spectrum near flux minimum. In general the bright 
phase spectra show stronger metal line absorptions than the 
corresponding continuum spectral class; Na I D at the stellar 

1 http://nedwww.ipac.caltech.edu/level5/Gray/Gray_contents.html 

2 http://www.sdss3.org/dr8/spectro/sspp.php 



TABLE 1 
J2339-5033 System Parameters 



Parameter 


ELC fit value 


RA (J2000) 


23:39:38.75 


DEC (J2000) 


-05:33:05.3 


P. (d) 


0.19309790±1 x 10~ 7 


T (MJD-TDB) 


55500.4833±0.0002 


Mi(M ) 


0.075±0.007 


M 2 (M Q ) 


1.40±0.04 


i(deg) 


57.4±0.5 


h 


0.90±0.01 


Ti (K) 


2800±50 


log[L x ](erg/s) 


33.5±0.1 



radial velocity is particularly strong at all phases away from 
inferior conjunction. 

Velocity shifts of individual exposures were computed by 
cross correlation against dwarf star spectral templates (from 
the SDSS spectral library) using the IRAF rvsao package. 
Templates from F0 through M5 were used. We found that K0 
templates provided the largest cross-correlation amplitude for 
most phases, despite continuum colors matching hotter spec- 
tral classes near maximum - we attribute this to the fact that 
the effective spectrum is composite. 

Near phase 0.5 the spectrum is increasingly dominated by 
molecular bands and the best cross-correlations (and effec- 
tive spectral class) move through M0 to M5. At these phases, 
the spectra are very red and improved correlation was often 
found when the data were trimmed by 500-1000A in the blue. 
Even at flux minimum we obtained cross-correlation coeffi- 
cients R > 3, except for a few poorly exposed spectra at large 
tracker offset (small effective telescope aperture). These were 
dropped from the radial velocity study. Typical correlations 
had i? ~ 15 — 25 and radial velocity errors 8-12 km/s . The ra- 
dial velocity points and errors are shown in the bottom panel 
of Figure 1, with symbol type indicating the highest correla- 
tion spectral class. 

3. BINARY SYSTEM PROPERTIES 

Normally for pulsar binaries we have a precision mass 
function and orbital ephemeris from the radio arrival times, 
which greatly restrict the range of viable model solu- 
tions. However, even without a pulse detection, we can 
still use our HET/SSO/WIYN photometry and spectroscopy 
to place preliminary constraints on the orbital parame- 
ters of the J2339— 0533 system, using the ELC code of 
lOrosz & Hauschildtl (12000b . This code includes the effect of 
companion illumination 'X-ray heating' and estimates light 
curves based on the low temperature 'NextGen' atmosphere 
tables. A table adapted for SDSS filters and white light pho- 
tometry was kindly generated by J. Orosz. As found by 
iRevnolds. et ail d200l for PSR B 1957+20, the low T eff at- 
mospheres were essential to allowing the model to produce 
reasonable light curves; the black-body approximation pro- 
duced models too bright at orbital minimum. For the fits we 
adjusted the companion unheated (backside) effective temper- 
ature T\ and Roche lobe fill factor /i, the pulsar heating flux 
L x , and the orbital semi-major axis a, mass ratio q, inclination 
i and phase of superior conjunction. 

The best fit parameters (Table 1, light curve and radial ve- 
locity models in Figure 1) indicates a cool companion, in- 
termediate inclination and moderate X-ray heating. This so- 
lution suggests a mass ratio q = 18.5 ± 1 and a relatively 
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massive companion Mi = 0.075 M®. The errors quoted are 
statistical only. Figure 1 also shows the companion center-of- 
mass (CoM) velocity; note the substantial decrease in ampli- 
tude and small shape distortions of the center-of-light (CoL) 
curve. We caution that, with a reduced \ 2 ~ 3, systematic 
errors are substantial and secondary fit minima are present. 

Overall, the good agreement with the 'black widow' pic- 
ture is encouraging. For all fits the models remain slightly 
too hot ( (/ too bright, i' too faint) at pulse minimum, suggest- 
ing amendments to the pass bands or the heating and albedo 
model are needed. Also the best fit epoch from the radial ve- 
locities alone shifts 8<j) ~ 0.005(5T = 0.001 d) from that of 
the full fit, implying that improvements in the spectral weight- 
ing could be helpful. Additional color photometry and espe- 
cially a pulsar mass function and phase measurement should 
tighten up the fits considerably. 



4. CONCLUSIONS AND IMPLICATIONS 

Our data also give some constraint on the distance and spin- 
down luminosity of the putative pulsar. At phase 0.25 the 
companion is half-illuminated, with T e f / » 4800 ± 200 K. 
Adopting the companion fit radius Ri = 0.24i? Q and scaling 
from a main sequence star gives an expected M r — 8.7 ± 0.3. 
The observed m r = 18.8 ± 0.2 then indicates a distance of 
1.1 ± 0.3kpc. The CXO flux gives a pulsar X-ray luminosity 
of w 3 x 10 31 erg s" 1 d^ pc . This is substantially less than the 

3 x 10 33 erg s _1 he ating effect ob served from the companion. 
Note, however that lBeckerl (120091) has found that pulsars have 
Lx ~ 10~ 3 i? so that we would infer a spin-down luminosity 
of - 3 x 10 34 d 2 ergs" 1 . Thus, again like PSR B1957+20, 
we find that a substantial fraction ~ 0.1 — 0.3 of the expected 
spindown luminosity goes into heating the companion. We 
can obtain an additional estimate of the pulsar spindown lu- 
minosity from the observed F 1 = 3 x 10~ n erg s _1 , since 

a heuristic 7-ray luminosity L 7t heu ~ ( E x 10 33 erg/s) 1 / 2 
has been found for other 7-ray pulsars (lAbdo et al.l 120101) . 
where L 7 — Att/qF^cI 2 and fa should be in the range 

0.7 - 1.3 Watt ers et all (|2009). This gives an estimate E w 
1.3 x 10 /j|d 4 ergs _1 . All of these estimates support the 

identification with a E ps 2 x 10 34 ergs _1 pulsar, somewhat 
less luminous and closer than PSR B 1957+20, but consistent 
with the other BW-type p ulsars recently discovered in the di- 
rection of Fermi sources (lRobertsll20nl 

We conclude that the 7-ray source is an energetic pulsar 
evaporating a low mass companion. Our preliminary pho- 
tometry and model fits suggest a somewhat larger companion 
mass than for the original black widow pulsar PSR B 1 957+20. 
Our present estimate of i ~ 57° implies that we should not 
see a 7-ray eclipse. However the substantial drop of the X- 



ray flux at phase 0.3-0.6 suggests a strong wind absorbing 
much of the X-ray emission. 

To date no 7-ray millisecond pulsar has been discovered in 
a blind search of the LAT photons, since the ~ms periods re- 
quire many trials in pulsar spin frequency and position space. 
Since most MSP are in binaries the additional trials needed to 
cover binary acceleration makes the searches prohibitively ex- 
pensive. One of our prime motivations in this project has been 
to reduce the binary phase space of trials needed to search 
J2339— 0533 and enable such a discovery. Our data tightly 
constrain the period and phase of the pulsar reflex motion and 
suggest an amplitude 02 sini rs 0.16 lt-s. These constraints are 
insufficient for a direct fold, but greatly restrict the number of 
binary trials needed. 

If a 7-ray (or radio) pulse detection is made, this system 
provides an excellent opportunity for study of the pulsar mass 
and binary evolution. The resulting precise mass ratio q 
will restrict the fit space, improving mass estimates for PSR 
J2339— 0533 and its companion. Further optical and X-ray 
spectroscopy of this bright system can measure the surface 
heating in more detail and probe the density of the evapora- 
tive wind. If highly eclipsed radio emission is discovered, 
this suggests radio beams closely coincident with the 7-ray 
emission. Conversely, if no pulsed radio emission is found 
while viewing at relatively clear orbital phases and folding 
at a known pulsar period, this could represent the first truly 
'radio-quiet' millisecond pulsar. As it is, we have built such a 
strong circumstantial case that it is fair to refer to J2339— 0533 
as a probable millisecond pulsar; all we are missing is the spin 
period itself. 
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